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We study, analytically and numerically, the classical and 
quantum properties of a nearly spherical 3D billiard. In par- 
ticular we show the appearence of quantum non ergodic be- 
haviour and of the deviations from Random Matrix Theory 
predictions which are due to the quantum suppression of clas- 
sically chaotic diffusion. 

PACS number: 05.45. -Kb 



In the ergodic theory of classical dynamical systems, 
billiards have played a fundamental role since they ex- 
hibit a clean and rich variety of dynamical properties, 
from completely integrable to true random motion. For 
this reason, since the early papers Q in which level spac- 
ing distribution was numerically studied, they became of 
primary importance for the theoretical and experimental 
analysis of the various quantum properties which are con- 
nected to different features of classically chaotic motion. 
An important question is to understand under what con- 
ditions the predictions of Random Matrix Theory (RMT) 
can be applied to a generic conservative quantum system 
In this context the quite surprising phenomenon of 
localization inside the energy shell has been shown to 
play a fundamental role This phenomenon implies 
quantum non-ergodic behaviour and deviations from the 
predictions of Random Matrix theory. Moreover a very 
interesting asymmetry appears between the structure of 
exact and unperturbed eigenfunctions. 

In a recent paper Q quantum localization was shown 
to take place in a stadium billiard and this result was 
confirmed also in a similar model The important 
question is now what happens in more dimensions. In 
particular we would like to understand whether or not 
this interesting non ergodic quantum behaviour can take 
place in a 3D billiard. 3D billiards are certainly more 
complicated than plane billiards. However, in some cases, 
they can be more close to experimental investigations |^- 
and to physical applications; we refer for example to 
dielectric resonators which may find applications in mi- 
crolasers and fibre-optic communications [ p^ 

In this paper we consider a wiggly 3D spherical billiard, 
namely a particle of unit mass and velocity v moving 
freely inside a closed 3D domain and bouncing elasti- 
cally off the boundary whose shape is given by a small. 



smooth but wiggly deformation of the unit sphere. The 
boundary of the billiard is described by the distance of 
the boundary point from the origin as a function of the 
unit direction vector n, rsin) = 1-1- ^f{n) {rp = 1). 
The small parameter e controls the size of the wiggles 
of the boundary while f{n) is an oscillatory function 
which is assumed to have unit average square modulus 
= (l/47r)/d2n[/(n)]2 = 1. The shape func- 
tion /(n) may be expanded in terms of spherical har- 
monics Yim{n) with a given maximal order L, f{n) = 
CimYim{n). The coefficients cim have been speci- 
fied by imposing three conditions: (i) decrease with 
increasing order i as c^^ oc l/t"^, (u) the shape of the 
billiard is invariant under the cubic symmetry group Oh, 
f{n) = f{Gn), G & Oh, (ni) the shape function /(n) is 
separable in cartesian coordinates, i.e. it can be written 



in a form f{n) 
cients a^. 



Ep=2 «p("x^ + + with coeffi- 



By calculating the distribution of maximal Liapunov ex- 
ponent, the classical dynamics of the billiard has been 
found to be almost completely chaotic for e > €c{L), 
where the critical chaos border ec{L) is a rapidly decreas- 
ing function of L. For the case L = 14 which will be 
taken in this paper, we found Cc = 0.001. For this value 
of the perturbation, most of phase space is covered by 
chaotic orbits. The billiard's dynamics can be described 
in terms of the angular momentum / — fx v and the unit 
direction vector ft — f/r of the forthcoming point of col- 
lision with the boundary. Labelling points of successive 
collisions by the discrete time variable j, one can write 
an approximate Poincare map for the dynamics in (I, ft) 
variables: 



{v^-l^)'^%xVf{n,)]+0{e'), (1) 



which is essentially four dimensional since we have two 



constraints, namely I 



J ■ ftj — and fi'j — 1. The compo- 
nents of Ij are the slow or adiabatic variables while fij 
is the fast variable (only one of its components is inde- 
pendent due to constraints). By averaging over the fast 
variable one can calculate the average drift of the angular 



momentum vector AL ^ I 
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From eq. (||) it is seen that the direction of angular mo- 
mentum vector undergoes a precession with a frequency 
vector field ^{l) = ©(eL^/^) which is an oscillatory func- 
tion of angular momentum on a scale vjL. 

Let us first analyze the behaviour of the modulus I of 
angular momentum. From we see that the average 
drift of I is zero (AZj) = (Z^+i - Ij) « (ijllj) ■ {^Ij) = 0. 
Then, above the chaos border e > ec, one typically ex- 
pects the diffusive behaviour {{Ij — ?o)^) ~ D{Iq)J where 
the diffusion coefficient can be estimated from (||) as 
D = 0{£^v'^ L'^). However, our numerical analysis reveals 
that the above diffusive behaviour takes place only up to 
time J — J* where J* « L. Instead, for J > J*, we find 
a slower, sublinear diffusion, which obeys the empirical 
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law {{I J - lo)') (X J" with < a < f [| 

For the variance of the component Iz of angular mo- 
mentum, both, the drift and diffusive terms are impor- 
tant ({Izj - Uo?) « {Izj ~ Izof + J{{{Alzo?) - (Al^of) 
w J^{n X lo)l + JD{lo). The first (drift) term 
is 0{J^e^v^L) and the second (diffusive) term is 
0{Je^v^L^), so the drift may be neglected w.r.t. dif- 
fusion for J < J* ^ L, while for J > J* the quadrati- 
cal drift term prevails [|l2|. It is not the purpose of the 
present paper to discuss the details of the classical dif- 
fusive mechanism. We would like only to stress that, as 
confirmed by our numerical analysis, the 'diffusion' of 
the direction of angular momentum is much faster than 
the diffusion of its modulus which is therefore the slowest 
variable in the system. 

The classical steady state equilibrium distribution Peq(0 
of the angular momentum can be obtained analytically by 
substituting (provided e > Cc) time averages with phase 
averages: 



I 

Pcq(0 - lim -Y,^tAl 

too 1 



n,0 



(3) 



Here Atj = 2(11^ — Ij^^^/v"^ are the time intervals be- 
tween collisions. As a measure of the width of the clas- 
sical distribution we introduce the quantity a{la,t) = 

{1/v) y^/o dlp{l, t)P - ( jp" dlp{l, t)lY with p{l,t = 0) = 

5{l — ^o)- For the equilibrium steady state we find 
(Teq = ct(/o,oo) = 0.2303. We have used espression (||) 
as an additional numerical test of ergodicity of our bil- 
liard for e > ec. 

Let us now turn to the quantum dynamics. We have 
solved the Schrodinger equation (V^ -I- fc2)4'fe(r) = 
with Dirichlet boundary conditions {fi — 1,E = k^/2). 
Eigenvalues /c^"^ and eigenfunctions ^'^cn) have been com- 
puted very efficiently using the scaling method invented 



by Vergini and Saraceno ||l3|l and already implemented 
for smooth 3D billiards in |14|]. Note that we have stud- 
ied only the states belonging to 1-dim fully symmetric 
irrep of the 48-fold cubic group Oh, so the Weyl for- 
mula which relates the wavenumber k with the sequen- 
tial quantum number J\f in this symmetry class reads (for 
small e) Af{k) « fc'^/(2167r). For sufficiently large k we 
expect that the statistical properties of eigenvalues and 
eigenfunctions of our billiard are well described by RMT. 
On the other hand, for sufficiently small k (or small e and 
fixed k) the perturbation theory is adequate and one can 
treat the problem as a small perturbation of a spherical 
billiard. We estimate the perturbative border ep(fc) as 
the critical value at which the change in the unperturbed 
levels, due to perturbation, is of the order of the spac- 
ing among levels, and therefore avoided crossings start to 
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appear. Note that the unperturbed eigenvalues fc, 
the sphere are (21 + l)-times degenerate with an average 
spacing (A/cq) ~ 8/A:. The derivative with respect to e 
of a level ko of a spherical billiard with an eigenfunction 
^^oir) can be expressed as an integral over the boundary 
of the billiard d,ko = -{l/2ko) J d^Sf{f/r)[dn'foir)? 
where 9„ is a normal derivative w.r.t. boundary. Then, 
using the explicit form for /(n) and \E'o(r), we can es- 
timate the level velocities deko « koy/g/{L + I) where 
g = 48 is the symmetry factor. Therefore we obtain: 



ep{k) 



(Afc) 8(L-fl) 



In fig. 1 we show the dependence of the eigenvalues on the 
perturbation parameter e. It is clearly seen that avoided 
crossings appear above the perturbative border which, 
for the parameters of fig. 1 (L = 14, fc « 180.0), is 
Ep « 0.0005. Above the perturbative border, namely for 
e > ep the quantum diffusion takes place according to the 
classical one. Our purpose is to try to understand if and 
under what conditions interference terms can suppress 
quantum diffusion, thus leading to localization and devi- 
ations from predictions of RMT. As it is known, quan- 
tum localization may take place after the Heisenberg time 
tn — 2'KdM / dE — fc/36 at which the quantum motion 
resolves the discreteness of the energy spectrum. More 
precisely, eigenstates are expected to be localized in the 
angular momentum coordinates if the Heisenberg time 
tu is less than the ergodic time tE which is the time it 
takes for the classical distribution to reach the equilib- 
rium state. Therefore, transition to ergodicity or delocal- 
ization takes place at tn = ts- If we assume, for the time 
being, that the classical motion is diffusive with diffusion 
coefficient D w e^L'^E, then the ergodic time is given 
by tE w {At)ll^^jD w 8/{3e^L^k) where (At) = 4/(3fc) 
is the average time between collisions. The condition 
tE = tn leads to the estimate for the delocalization bor- 
der ei: 
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However, due to anomalous difFusion, the ergodic time 
will be longer than the one given by the above esti- 
mate and this will lead to a delocalization border even 
larger than (^). For example, for the data of figs. 1-4 
(L = 14, /c « 180) the numerically estimated localization 
border is e/ = 0.010. 

Billiard's cigcnfunctions ^'fe(r) can be expanded in 
terms of eigenfunctions of a perfectly spherical bil- 
liard {r\nml) = {V2/ji+i{^in))jii^inr)Yi^{n), l^-fe) = 
^j^j^(nZm|^'fc)|nZ77i). We define the angular momentum 
distribution of an eigenstate '^k with an eigenvalue k as 
^^kil) = p ^(;) J2 \{nlm\'i'k)\^- Since we have divided by 

the classical stationary distribution Pcqil) then, for com- 
pletely delocahzed, ergodic states, hk{l) should approach 
a constant (apart from fluctuations). In fig. 2 we show 
the distributions hk{l) for three typical eigenstates of the 
billiard at three different values of parameter e. Eigen- 
states (a) and (b) are below the delocalization border 
and are exponentially localized while (c) is an ergodic 
extended eigenstate. However, eigenstates {nl'm,\'i>) are 
extended in the Iz = rn variable as we expect as a con- 
sequence of the much faster classical diffusion of 1^ w.r.t. 
I. 

For a comparison with theoretical predictions, in anal- 
ogy with the classical case, we measure the size of 
eigenfunctions by the quantity (7^ — jii"^ k\l'^\'^ k) ~ 
(^ffcl^l^fe)^)^/^. We expect that for ergodic (extended) 
eigenstates, at should agree with the width of the classi- 
cal steady state angular momentum distribution acq = 
0.2303. Instead, for localized states, we expect that 
(cTfc) ~ {a{lo,T))ig where r is of the order of Heisen- 
berg time. In order to suppress fluctuations the quan- 
tum width is averaged over a sufficiently large number of 
consecutive eigenstates while the classical width is micro- 
canonically averaged over the initial angular momentum 
Iq. In order to take into account that, as is empirically 
known, the width of eigenstates is approximately one half 
of the width of steady state we choose t = <///2 = fc/72. 
The comparison between (ak) and {a{l,tH /'2,))i is shown 
in fig. 3 for a fixed energy window. The agreement is sur- 
prisingly good. Moreover we have numerically checked 
that the above relation also holds for different energy 
windows with k ranging between k w 120 and k ~ 260 
and for fixed e = 0.003. 

We have also calculated in the different regimes the 
most commonly studied statistical spectral properties, 
namely the nearest neigbour level spacings distribution 
P{S) and the number variance T,'^{E). The distribu- 
tion P{S) has been found to be characterized by the 
power-law level repulsion P{S) oc , S ^ 1, where /3 
smoothly increases from zero to one as a function of e or 
k on moving from the perturbative to the ergodic regime. 
The functional relation (3{{(Jk)) is still to be systemati- 
cally investigated. However, in the perturbative regime 
e < ep{k), f] Ri but the tails of P{S) decrease much 
slower than Poissonian cxp(— S*) due to strong degenera- 
cies at e = 0. Moreover the convergence towards GOE 



distribution with increasing e or fc is found to be very 
slow and in agreement with previous results on differ- 
ent models [3]. In this respect, the S^(£^) describing the 
long range correlations has been found to be more sen- 
sitive. In fig. 4 we plot S^(£') computed in the window 
170 < fc < 180 and for e = 0.003, 0.010, 0.016. 

In this paper we have shown, on a model of classically 
chaotic 3D billiard, the appearence of the quantum dy- 
namical localization phenomenon which results in the ap- 
proach of the quantum distribution to a localized, non 
ergodic, steady state. This phenomenon, though in a 
different context, was discovered long ago [|l^ and it is 
related to the existence of different time scales of classical 
and quantum motion ||l6|l . 

Finally we would like to remark that the 3D nearly 
spherical billiard studied here can be a realistic model 
for a new class of optical resonators | [lO| |: in particular 
the localization in angular momentum can lead to an 
increase in the photons lifetime. 
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FIG. 2. Angular momentum probability distribution hk{l) 
of three eigenstates with k « 180 and different e: (a) strongly 
exponentially localized cigenstatc for e = 0.001 with eigen- 
value k = 180.3009, (b) localized cigenstatc for e — 0.003 
with eigenvalue k = 179.8013, and (c) extended (crgodic) 
eigenstate for e = 0.016 with eigenvalue k = 179.8611. Note 
that the probability distribution of the state (a) is divided by 
a factor 10® and probability for the state (b) by a factor 10^ 
w.r.t. state (c). 
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FIG. 1. Dependence of quantum eigenvalues on the per- 
turbation parameter e in the window 179.68 < k < 180.32. 
This interval contains approximately 90 levels. 
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FIG. 3. Comparision of the average width (<7k) of angular 
momentum distribution of quantum eigenstates (bullets) and 
the corresponding widths {a{lo,tH of classical distribu- 
tions (diamonds) after half of Heisenberg time for different 
values of the perturbation parameter e and fixed L = 14. The 
quantum distributions have been averaged over 90 consecutive 
states around k w 180 
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FIG. 4. Number variance statistics S^(£) for three spec- 
tral stretches in the window 170 < fc < 180 corresponding 
to three different values of the perturbation parameter (a) 
e = 0.003, (b) e = 0.010, (c) e = 0.016. Each window con- 
tains about 1500 levels. The dashed curve is the logarithmic 
GOE number variance while the dotted curve is the linear 
(Poissonian) S^(i5) statistics. The relevant energy scales for 
this plot arc: mean level spacing AE = 72n/k ~ 1.26; Thou- 
less energy Ec = l/ts ~ S ior curve (c); the bouncing energy 
Eb = 2/{t) m 270 where {t) is the average bouncing time. 
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